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Properly modulated Hatband lattices have a divergent density of states at the Hatband energy. 
Quasiperiodic modulations are known to host a metal insulator transition already in one space 
dimension. Their embedding into Hatband geometries consequently allows for a precise engineering 
and Hne tuning of mobility edges. We obtain analytic expressions for singular mobility edges for two 
Hatband lattice examples. In particular, we engineer cases with arbitrarily small energy separations 
of mobility edge, zeroes, and divergencies. 
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I. Introduction 


The phenomenon of wave localization has been in¬ 
tensively studied since its prediction in 1958 [T], where 
complete localization was proved in the case of a one¬ 
dimensional (1-D) chain defined over a random poten¬ 
tial. Moreover, it was shown that the 3-D case allows for 
an energy-dependent transition from localized to delo¬ 
calized eigenstates. The transition has been since coined 
metal-insulator transition (MIT). The critical energy Ec 
is called a mobility edge; in general, it depends on and 
varies upon changes of the control parameters of the given 
model [2]. Interestingly, an MIT can also be realized in 
one-dimensional settings with sufficiently correlated dis¬ 
order potentials [3]. 

In 1980, Aubry and Andre proved the existence of the 
MIT for a I-D chain defined over a specific quasiperi¬ 
odic potential [4]. This MIT occurs at a critical value 
(Ac = 2) of the onsite potential’s strength A, and sep¬ 
arates the metallic phase A g]0,2[ from an insulating 
phase A g]2,+oo[. This remarkable result was fully un¬ 
derstood via the principle of duality, in which a partic¬ 
ular Fourier transformation relating eigenmodes and en¬ 
ergy spectra allows for a direct functional equivalency be¬ 
tween momentum space and its transform counterpart. 
This equivalency is energy-independent: upon crossing 
the critical value Ac all eigenstates turn from localized 
to extended, regardless of their eigenenergy. The ap¬ 
pearance of a mobility edge is thus avoided. Analytic 
results have been discovered in the last decade regard¬ 
ing the topological Cantor structure of the spectrum [5] 
and its Lebesgue measure [6]. Furthermore, for each dif¬ 
ferent regime (insulating, metallic and critical) different 
spectral decompositions have been found [THS] . Model 
generalizations were reported, e.g. quasiperiodic systems 
constantly maintained at criticality nnmi], bichromatic 
quasiperiodic lattices displaying mobility edges [121 US] , 
and completely localized quasiperiodic models [14]. Cor¬ 
related metallic states have also been observed in the in¬ 
sulating regime, for the case of two interacting particles 
within a 1-D Aubry-Andre chain [T5|. In another recent 
work [16] , a suitably modified quasiperiodic potential was 


shown to produce a mobility edge expressable in an ana¬ 
lytic form - a property which we will take to new limits 
using Hatband topologies. 

Wave propagation on lattices with Hatband topologies 
is characterized by the existence of horizontal (fiat) bands 
in their band structure. Known in condensed matter, 
this model class has gained great interest in the scientific 
community, due in part to experimental realizations in 
optical lattices and paraxially approximate light propa¬ 
gations [Mini- Recent theoretical discoveries have also 
considered the presence of a disordered potential HEI] 
and nonlinearity [22] . An innovative procedure detangles 
Hatband states from dispersive ones [23] . This allows one 
to inspect specific features of the models as they relate 
to the choice of the onsite perturbations, and also sug¬ 
gests specific potential correlations. In the present work, 
this detangling technique of local rotations [23] is applied 
as an extension of [24]; in particular, regarding the pre¬ 
liminary finding of a MIT occurring in a Hatband lattice 
under quasiperiodic Aubry-Andre perturbation. 

The present paper has the following structure: in Secjn] 
the general features of Hatband topologies are introduced 
that define two particular models (cross-stitch and dia¬ 
mond lattices), a quasiperiodic Aubry-Andre onsite per¬ 
turbation is defined, and the coordinate transformation 
that allows rotation into Fano defect lattices [23]. In 
Sec |III||lV| our findings for the cross-stitch and diamond 
lattices are respectively presented: for both, two distinct 
chain correlations are discussed. Where applicable, the 
exact mathematical expression obtained for the mobility 
edge is jointly shown with numerically obtained transi¬ 
tions for these particular onsite correlations. 


II. Flatband Topologies 

Consider the eigenvalue problem of a generalized tight- 
binding model 

Elpn = eni’n “ V'lpn “ 7(V’n-l + V’n+l) • ( 1 ) 

For all n G Z, each component of the vector = 
( 7 /;^,..., 7 /^^)^ represents a site of a periodic lattice, while 
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the set of sites represented by ipn is the n-th unit cell. 
The real matrix V defines the geometry of the unit cell, 
while the real matrix T describes hopping to neighboring 
cells. At each of the lattice’s Tth leg {'iph}n^ an onsite 
perturbation is defined. The unit cell perturbation 

€n of Eq. Q is thus given by the diagonal square matrix 
e„ = diag(e“,e^,...,e^). 

The model geometry is contained in the matrices T and 
E, which are then used to derive the dispersion spectrum 
via the Bloch solution on an unperturbed 

crystal = 0. Platband topologies are models in which 
this crystalline case exhibits at least one band indepen¬ 
dent of k - such a band is dispersionless, or “fiat”. Eigen- 
modes corresponding to this fiatband energy are (usually) 
compact localized states (CLS), i.e. modes whose ampli¬ 
tude is nonzero only across a finite number of sites [23] . 
The fiatband topology class, t/, is then defined as the 
minimum number of unit cells the CLS occupies [23] . 

In this paper we consider two lattice topologies - the 
cross-stitch and diamond lattices. The former, shown in 
the upper left of FigHJ is defined for a unit cell pJn = 
{an^ bn)^ with a 2 X 2 peturbation matrix e^. This yields 
for Eq. 0 the following matrices 

'), Tcs=(;;). ( 2 ) 




FIG. 1. Left: The cross-stitch lattice - the grey shaded region 
indicates the unit cell. Right: The transformed Fano defect 
lattice of Eq. ^ . 


Likewise for the diamond lattice, as shown in the upper 
left of Eigj^ the unit cell is pJn = {ctn^bn, Cn)'^ with a 
3x3 perturbation matrix e^. In this case, the matrices 
in Eq. 0 are 

/0tl\ /ooo\ 

Vdc = { t 0 1 \ , = 0 0 0 . (3) 

Vi 1 oy Vi 1 V 



FIG. 2. Left: The diamond lattice - the grey shaded region 
indicates the unit cell. Right: The transformed Fano defect 
lattice of Eq. (20). 


In the unperturbed crystal = 0, the dispersive bands 
are 


{ —t — 4 cos k, Cross-Stitch, 

— - (t ± \/t^ 16 cos k + 16) , Diamond. 

Additionally, both models contain a fiat band at E{k) = 

t. 

Associated with the fiatband energy, a resulting CLS 
can be constructed: pJn = (I 5 ~l)^^n,no/v^ (cross- 
stitch) and pJn = (I 5 “I, 0 )^^n,no /\/2 (diamond). Note 
that both CLS are contained within a single unit cell. 
Therefore, according to the definition previously stated, 
both lattices are fiatband models class U = 1. Coordi¬ 
nate transformations local to the unit cells rotate these 
lattices into a Eano defect form [23]. Eor the cross-stitch, 
the rotation is defined by the real matrix Ucs 

(^P£^=Ucs^Pn, Ucs = ^l^l\y (4) 

Similarly for the diamond lattice the transformation is 
defined by the real matrix Udc 

(pA ^ ^ 1 C ^ M 

/n = UdC V’n, ^ 1 -1 0 . (5) 

\Cn) 0 V 2 ) 

Lastly, this local coordinate transformation must also ro¬ 
tate the onsite perturbation. Eor both lattices, this gives 

e± = (e“±e^)/2. (6) 


The effect of quasiperiodic Aubry-Andre perturbations 
on these two topologies is the focus of the present work. 
Eor both lattices, the onsite perturbations are de¬ 
fined as independent Aubry-Andre potentials 

e\ = \iCos[2TT {an^ Oi)] , (7) 

for the i = aA (cross-stitch) and i = a, 6 , c (diamond) 
legs. The parameters are positive real values control¬ 
ling the perturbative strength, Oi is the phase-shift, and 
a is an irrational number (here set to the golden ratio) 
called the ineommensurate parameter. Without loss of 
generality, the a-leg phase can be zeroed {Oa = 0). We 
also set the leg potential strengths equal to each other 
Ai = A. 

Erom Eq. 0 , notable correlations between the a-leg 
and 6 -leg perturbations appear and will be object of our 
studies for both models; namely 

Symmetric; e;; = 0 
Antisymmetric: = 0 . 

Since the a-leg phase has been zeroed, from Eq. 0 these 
two correlations are obtained solely from the 6 -leg phase, 
e.g. 6^5 = 0.5 ( 6^5 = 0) for the antisymmetric (symmetric) 
case. We start the analysis of these models with the 
cross-stitch in Sec |III[ and then with the diamond lattice 
in Sec lIVl 
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III. Cross-Stitch Lattice 


By Eqs.(|4|6|), the cross-stitch lattice transforms into 


{E ^t)pn= PnEe^ fn- 2 {pn-1 + Pn+l) , 


{E -t)fn= e+ fn + e„ . 


( 9 ) 


This results in a Fano chain, as shown in the right of 
FiglU The local rotation yields a dispersive coordinate 
Pn and a compact Fano coordinate /^. The sequence 
describes onsite perturbations of both Pn and /n, while 
the sequence e~ couples the dispersive to the Fano coor¬ 
dinate within the rotated unit cell m- Solving for the 
Fano coordinates in the second equation above, we 
obtain a new equation for the dispersive portion 


{E -\-t)pn — 


-^2 


(en) 


(E-t)- ei 


Pn -2{pn-l +Pn+l) ■ 

( 10 ) 

The reduced topology assumes the tight-binding form. If 
eigenmodes are exponentially localized, their asymptotic 
decay is e~ ^. The rate ^~^{E) is the inverse local¬ 

ization length of a localized state at eigenenergy F G M, 
found by applying the recursive iteration 


rHE,X)= lim ,, 

^ ^ M->+oo M 


M 

\/f ^ 


n=l 


Pn-\-l 


Pn 


( 11 ) 


for any given potential strength A. We will use this 
method in all the numerical computations of the two 
models’ localization lengths, for M = 10^. The energy 
E in Eq.(lO) will be numerically found from an exact di- 
agonalization of a finite lattice of N = 512 unit cells. In 
all the figures of the paper, if the recursive iteration con¬ 
verges to a finite value (chosen [25] here as ^ < A^/10), the 
datapoint (F, A) is declared a localized state and plotted 
in blue. Otherwise if the iteration diverges, the datapoint 
is declared an extended state and plotted in red. 

A. Symmetric Case: Metal-Insulator Transition 

We analyze first the symmetric case e“ = 0, obtained 
for Ob = 0.0. Eq. reads 


{E -\-t)pn — Pn — 2 {Pn-1 + Pn+l) 5 
{E -t)fn= e+ fn 


( 12 ) 


with The two sets of states Pn and fn decou¬ 

ple and generate two independent spectra, respectively 
labeled ap and a/. The parameter t then simply oper¬ 
ates as a shift parameter, translating ap and af relative 
to each other by 2t. 

The dispersive states Pn are described by an Aubry- 
Andre chain, displaying a MIT at Ac = 4. The af states 
keep their compact feature, but the degeneracy of their 
eigenenergies is now removed - these eigenenergies are 
given by F = e+ In Figj^ we plot the spectrum from 
Eq.( 12 ), as a function of A. In this symmetric case, the 


spectra ap and af are independent, since the Pn and the 
fn coordinates decouple. 

For every potential strength A > 0, the Fano states 
spectrum af = is equidistributed within the in¬ 

terval [t — A,t + A]. In FigJ^we indicate its boundaries 
by dashed lines. For the dispersive spectrum cr^, the lo¬ 
calization length is numerically found with the recursive 
iteration 0 , and the localized phase (red) is demar¬ 
cated from the extended one (blue). At the critical val- 
uce Ac = 4, all dispersive states switch from extended to 
localized. In this case, there is no mobility edge. 



FIG. 3. Symmetric Case: The dispersive spectrum {ap) of 
the cross-stitch lattice, for e“ = 0 and t = 0. The Fano 
state spectrum af is omitted, but its boundaries indicated 
by black dashed lines. The black line represents the MIT at 
A = 4, which clearly separates extended states (blue) from 
those localized (red). 


B. Asymmetric Case / 0: Numerical Evidence 
for Mobility Edge 

Breaking the symmetry e“ 7 ^ 0 {Ob 7 ^ 0) of the Fano 
chain Eq. effectively couples the dispersive states Pn to 
their compact Fano counterparts fn- Therefore, the self¬ 
duality is lost, and the two independent spectra {apj) are 
now joint. Nevertheless, we expect a transition between 
localized and extended states via an energy-dependent 
mobility edge. In Fig we plot the spectrum of the lattice 
in the asymmetric case for Ob = 0.25 and t = 0. A 
mobility edge is clearly observed separating the localized 
regime (red) from that which is extended (blue). 

C. Antisymmetric Case {ef = 0): Analytic Mobility 

Edge 

Among all the possible non-symmetric cases obtained 
for ^5 7 ^ 0, the antisymmetric one Ob = 0.5 deserves spe- 
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FIG. 4. Asymmetric Case: Spectrum of the cross-stitch lat¬ 
tice, for 0 — 0.25 and t = 0. The extended (localized) portion 
of the spectrum is shown in blue (red), while the color bound¬ 
ary is a mobility edge approximation of the spectrum. 


cial attention. In this situation = 0 and Eq.(|^ trans¬ 
forms into 


(-E “h t) Pji — ‘lijPfi—x + 5 

{E -t)fn = e~Pn. 


It follows that all Hatband states are expelled from their 
unperturbed energy position E = t. Since e~ ^ 0, from 
the second equation of (13) it follows that at the Hatband 
energy 


Pn = 0 


Pn — 0 : 


(14) 


Then, from the Hrst equation of (13) we conclude 

fn = 0 => /n = 0. 


(15) 


Only the trivial state {pn, fn) = (0,0) satisifies Eq. (HI 
exactly at the Hatband energy E = t. In Fig|^we plot 
the spectrum for this antisymmetric case for t = 0. We 
again observe a mobility edge. However, we can now even 
identify its analytical form, as observed in [24]. Indeed 
the dispersive part, Eq.([TQ|), reads 

(e-)2 

{E ^t)pn = Pn - 2(Pn-l + Pn+l) • (16) 

— t 

By trigonometric bisection 

(e“)^ = cos^(27r(an) ~ ^ + cos(47r(an)] . (17) 



FIG. 5. Antisymmetric Case: Spectrum of the cross-stitch 
lattice, for = 0 and t = 0. The extended (localized) spec¬ 
tral portion is shown in blue (red). The color boundary is 
an approximation to the spectral mobility edge, however the 
black line is the analytical form of Eq.([T^ - a good agreement 
is observed. 


Substituting Eq.([T7|) back into the previous equation, 
we obtain 


-E Pn — 


A^ 


where E := 


4{E -1) 
E 1 


cos(47r(an) — {pn 
A^ 

~ 4{E-ty 


1 +Pn+l), 


(18) 


Therefore, the model becomes an Aubry-Andre chain 
eigen-equation, but with onsite perturbation strength de¬ 
pending on A and E. From |4|, the MIT occurs when 
the potential strength is twice larger than the hopping 
strength. Imposing that condition, an analytic expres¬ 
sion is found for the mobility edge, Ac(Ec): 


4(Ee - 1) 


= 2 


Ae(Ec) = • (19) 


Note that for Ec = the critical potential strength Ac 
vanishes in a square root manner, where the previously 
discussed lack of states at the Hatband energy E = t oc¬ 
curs. The analytic curve of the mobility edge is plotted in 
Figj^ displaying excellent agreement with the numerical 
result. 
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IV. Diamond Lattice 


Under Eqs.([^[^, the diamond lattice’s Eq. 0 become 

(-S' “1“ Pn ~ ^nPn fn “1“ Cn+l)^ 

(-^ ~ t) fn — fn + Pni (^ 0 ) 

{E - el) Cn = -V2(pn-1 Pn), 

as illustrated graphically in the right plot of Eig|^ 
Expressing the fn and Cn variables through the Pn ones 
we reduce these equations to a tight-binding form which 
contains the Pn coordinates only: 


{E^t)pn — 


et-\- 




{E-t)e 


: + 


+ - 


E — E — e! 


n+l 


E 


rPn-l + 


E 


Pn+l • 


PnE 

( 21 ) 


■^n+l 


We use this tight-binding form to numerically obtain the 
eigenvalue spectrum and consequently the localization 
length using the iterative recursion Eq.(pd]). In EigJ^ 
we plot the diamond lattice spectrum and the mobility 
edge for ^5 = 0.5, Oc = 0.05 and t = 0. 



FIG. 6. Spectrum of the diamond lattice, for Ob = 0.5, Oc = 
0.05 and t = 0. The extended (localized) spectral portion 
is shown in blue (red), while the color boundary is again a 
spectral mobility edge approximation. 


We notice that in Eq.( 21 ) the coefficients of the hop¬ 
ping terms depend on the c-chain onsite perturbation e^. 
Therefore, in general all hopping terms are perturbation- 
dependent. In order to arrive at an Aubry-Andre form 
of Eq.(21) we set constant onsite energies on all c sites: 

= iF G M. Then an extended state V exists at en¬ 
ergy E = K, regardless of the other control parameters 


in Eq.(20): 


E = K , Cn = {-1)^ , fn=Pn = 0. (22) 

The state’s amplitudes reside only on the c sites (see 
Eigj^. The existence of this extended state is not af- 



FIG. 7. Extended state V at energy E = K on the diamond 
lattice in case of constant onsite potential = A on the 
c-chain (up to normalization). 


fected by the perturbation strength A, the hatband en¬ 
ergy A = t, and any specihc correlation of the on¬ 
site potential. Therefore if the model admits a mo¬ 
bility edge curve Ac (Ac), it follows that it will diverge 
Ac (Ac = A) = oo, yielding a singularity. 

A. Symmetric Case: Analytic Mobility Edge 


We consider hrst the symmetric case = 0, obtained 
for Ob = 0.0. In this situation, Eq.([^ reads 

(A + t) Pn = ^n Pn “I y^i^n + ) , 

iE-t)U =4fn, (23) 

{E - K)Cn = -V2{pn-1 + Pn) ■ 

The fn states decouple from bothp^ and Cn states, gener¬ 
ating two independent spectra a/ and cr^^c- The hatband 
energy t shifts the two energy spectra of 2t relative to 
each other. Substituting the c variables by the p ones we 
obtain the Aubry-Andre form 


~ (A - A)A 

Epn = - - - COs(27r(an)pn E Pn-l +Pn+1, 

, ~ (A + t)(A-A) ^ 

where A :=-2. 


(24) 


Note the onsite potential strength is now dependent on 
A and A. Imposing the equality between the potential 
strength and the Aubry-Andre critical value, we arrive 
at the mobility edge 


{Ec - K) 


= 2 


Ac(^c) 


4 


Ac-A 


(25) 


The mobility edge curve diverges at A = A due to the 
existence of the delocalized state V. We plot this mobility 
edge curve in Eigj^and observe very good agreement with 
the numerical results. 


B. Antisymmetric Case: Analytic Mobility Edge 

We consider the antisymmetric case e+ = 0 obtained 
with Ob = 0.5, and with = A. Eq.([^ reads 

{E + t) Pn — ^n fn E^i^n E Cn-\-l) , 

(A -t) fn = e~Pn , (26) 

{E - K)Cn = -\/2{pn-i + Pn) ■ 
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FIG. 8. Symmetric Case: Dispersive spectrum (cTp^c) of the 
diamond lattice, for e^=0, t = The Fano state 

spectrum (cr/) is omitted, but its boundaries are indicated 
by black dashed lines. The extended (localized) states of the 
dispersive spectrum are shown in blue (red) - the bound¬ 
ary between these is a mobility edge approximation, in good 
agreement with its analytical form of Eq. (25), shown as a solid 
black line. Note that at E = K, the mobility edge diverges 
to infinity. 



FIG. 9. Antisymmetric Case: Spectrum of the diamond lat¬ 
tice, for Cn =0, K = = —1. The extended (localized) 

spectral portion is shown in blue (red). The boundary be¬ 
tween is an approximation of the spectrum’s mobility edge - 
in good agreement with the analytical form of Eq. (28), plot¬ 
ted as a black line. Not the mobility edge curve zeroes at 
E = t, as well as diverges at E — K. 


For t 7 ^ iF, all eigenenergies are expelled from the flat- 
band energy E = t. In Figj^we plot the spectrum for 
this antisymmetric case for t = — 1 and K = 1. We derive 
an analytical expression of the mobility edge by reducing 
Eq.(26) to an Aubry-Andre form for the Pn coordinates: 


Epn = 


where E := 


E-K 


2 2{E-t) 

E-K 


{E + t) 


cos(47rQ;n) + p„_i + p„+i, 
A2 


2{E-t)_ 


2 . 


(27) 


The condition for the MIT yields 


E-K XI 

2 2{E-t) 


2 Xc{E) = 2 



E-t 

V 

E-K 


(28) 


In Eq.(28), the mobility edge curve Xc{Ec) diverges to 
infinity at E = K, in correspondence to the delocalized 
state V. The curve also displays a zero at E = which 
corresponds to the lack of any states at the hatband en¬ 
ergy [26]. The mobility edge curve of Eq.(28) is plotted 
in Figand show excellent agreement with the numerical 
data. 


V. Conclusion 

Flatband topologies are characterized by macroscopic 
degeneracy at the flatband energy. General perturba¬ 
tions of these topologies lead to a removal of the degen¬ 
eracy, yet keeps a high density of states and a bunch¬ 
ing of the renormalized and hybridized states around the 
original flatband energy. This has especially dramatic 
consequencies for quasiperiodic Aubry-Andre form per¬ 
turbations. The flatband energy now hosts a zero of 
a mobility edge curve Xc{Ec). When approaching this 
zero the density of states grows, and the spatial extent 
of the eigenstates drops, making them more localized. 
For specific symmetries of the applied quasiperiodic po¬ 
tential, the dependence Xc{Ec) is obtained analytically, 
confirming the predicted zero, and further proving the 
strict nonexistence of any state at the former flatband 
energy. Some flatband topologies allow the existence 
of completely delocalized eigenstates at certain energies 
E = K. This leads to even more complex mobility edge 
curves which allow for a coexistence of zeroes and diver¬ 
gencies of Xc{Ec). Possible future topics of study include 
extension to U > 1 topological classes and higher dimen¬ 
sional flatband models. It is our hope that the use of flat- 
band topologies contributes interest to tunable mobility 
edges, e.g. by those realized in graphene EH, monolay¬ 
ered dichalcogenides [28], or vandium dioxide films [29] . 
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